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Abstract 

Recently, there has been a significant interest in developing cooperative sensing systems for certain 
types of wireless applications. In such systems, a group of sensing nodes periodically collect measure- 
ments about the signals being observed in the given geographical region and transmit these measurements 
to a central node, which in turn processes this information to recover the signals. For example, in 
cognitive radio networks, the signals of interest are those generated by the primary transmitters and the 
sensing nodes are the secondary users. In such networks, it is critically important to be able to reliably 
determine the presence or absence of primary transmitters in order to avoid causing interference. The 
standard approach to transmit these measurements from sensor the nodes to the fusion center has been 
to use orthogonal channels. Such an approach quickly places a burden on the control-channel-capacity 
of the network that would scale linearly in the number of cooperating sensing nodes. In this paper, 
we show that as long as one condition is satisfied: the dynamics of the observed signals are sparse, 
i.e., the observed signals do not change their values very rapidly in relation to the time-scale at which 
the measurements are collected, we can significantly reduce the control bandwidth of the system while 
achieving the full (linear) bandwidth performance. 



*A preliminary version with a subset of the results was presented at the Allerton Conference on Communication, Control and 
Computing, Monticello, IL, Sep. 2010. 
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I. Introduction 

Cooperative sensing is a promising technique that has received a lot of attention recently due 
to necessity for reliable decision-making. Here, a group of wireless nodes collects measurements 
about the some signals of interest. The observations are then reported to a central node, which in 
turn applies some appropriately-chosen decision rule to recover the signals. Such an architecture 
finds application in two closely-related and well-studied settings. 

The first setting is (de-)centralized detection in sensor networks - a group of sensor nodes 
commanded by a fusion center - where the observed signals could be temperature or humidity 
readings, military targets, etc. Much work (see flU — [SJ and references therein), both theoretical 
and practical, has studied and optimized various aspects of the cooperative sensing architecture, 
e.g., quantization functions at the sensing nodes, the combining process and decision metric at 
the fusion center, the communication scheme between the sensing nodes and the fusion center, 
etc. This is a rich area of research and hence, our references are far from a comprehensive list. 
The second setting is a cognitive radio network where the observed signals are other unlicensed 
(secondary) or licensed (primary) transmitters. While the ideas in this paper are applicable to 
both the aforementioned settings, we review some literature from the latter as it is most relevant 
to the technical developments proposed here. 

A. Spectrum sensing in cognitive radio networks 

Cognitive radio technology, introduced by J. Mitola 0, is a promising solution to the growing 
scarcity of wireless spectrum, one that can potentially increase the spectrum utilization efficiency 
as recognized by the FCC [0, 0. Traditionally, a portion of spectrum is allocated or licensed 
for use by a specific group of users by regulatory agencies. On the other hand, cognitive radio 
networks call for cognitive (unlicensed/secondary) users to operate on the same frequency band 
as the primary licensed users while attempting to access the spectrum seamlessly. In other words, 
the cognitive users need to adjust their operating parameters to guarantee minimal impact to the 
primary licensed users. For example, DARPA's Next Generation program [JSj has been interested 
in developing spectrum sensing techniques that prevent interference to existing occupants of the 
frequency band. To increase spectral efficiency, the FCC recently opened up TV whitespaces 
(54 MHz - 806 MHz) for unlicensed use JH. While the ruling calls for access to a central 
database to determine TV band availability, the traditional use-case for spectrum sensing, we 
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believe it still plays an important role in providing acceptable quality-of-service (QoS) over 
unlicensed bands. For example, a database provider could provide multiple classes of service 
to a querying secondary user where the higher class would be provided with a per- TV-band 
estimate of the interference statistics from other unlicensed users in the band. Spectrum sensing 
might be necessary in this context in order to collect such statistics. 

Cooperative or collaborative sensing in the context of cognitive radio networks consists of a 
group of secondary users or specially-placed sensor nodes that collect measurements over some 
sensing time window about the activity of the primary user and transmit these measurements 
to a central node. The central node or fusion center may often be the cognitive base station. 
Cooperative sensing techniques can be classified appropriately based on the type of decision 
metric used at the fusion center. We focus on one of the simplest schemes that relies on energy 
detection. Most prior work on energy-detection-based cooperative sensing (see lfT0ll -[fT4 l| and 
references therein), the signal is typically a single binary hypothesis modeling a system with 
one primary user that is either ON or OFF. In cooperative sensing, the measurements from each 
cognitive sensor node are typically linearly weighted and combined to form a decision statistic 
that is in turn compared against a threshold to produce a binary output. The paper by Ganesan 
and Li ifTOl constitutes one of the earliest contributions that establish the gains in detection 
(of a single primary user) probability due to cooperation in the presence of channel fading. 
Peh et al. |fTT| and Lee 031 study the tradeoff between sensing time and throughput of the 
cognitive network since longer sensing times lead to higher primary detection probabilities but 
lower secondary throughputs. Quan et al. Ifl2ll further optimize the linear weighted combiner 
studied by Ganesan and Li ifTOll by choosing weights that maximize the detection probability. 
While the above authors designed soft-combining systems where the sensors typically report their 
measurements to the fusion center without any processing, Li et al. lfT3l consider hard-combining 
where local decisions are made at each sensor, and design a voting rule to produce the final 
binary decision. Finally, Mishra et al. lTT5l compare the performance of hard-/soft-combining 
and study the effects of correlated fading on the detection probability. We refer the reader to 
a recent survey papers by Akyildiz et al. Ifl6ll and Yucek lUTl for a more comprehensive list 
of references along with further discussion on the other types of cooperative spectrum sensing 
such as cyclo stationary feature detection. 
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B. Our contributions 

In almost all of the reviewed literature, the standard approach to transmit the measurements 
from the sensor nodes to the fusion center has been to employ orthogonal channels. This means 
that the control bandwidth demanded by the standard cooperative sensing scheme scales linearly 
in the number of sensing nodes. This immediately places a significant burden on the bandwidth 
requirements of the network and might be impractical in some scenarios as recognized by 
Akyildiz et al. lfT6ll . In this work, we focus reducing the amount of control bandwidth required 
by the system by considering the structure that may often exist in the dynamics of the observed 
signals. Thus, we immediately build on and extend past models by adding a temporal dimension 
to the observed signal and by considering multiple discrete signals instead of a single binary 
hypothesis. 

The key idea we exploit in this paper is the following: while the observed signals at a snapshot 
in time in general lives in some arbitrary high dimension, often, there is a time-scale separation 
between the sensing time window and the behavior of the observed process. For example, TV 
transmitters would turn ON/OFF on a significantly slower time-scale (in the order of minutes) 
than sensor measurement windows (in the order of milliseconds). Under these circumstances, the 
dynamics of the observed signal vector is likely to be sparse. It has long been known, and recently 
popularized under the name of compressive sampling, that whereas N linear measurements 
are required to reconstruct a vector or signal in M. N , if it is S-sparse (i.e., it has S non-zero 
coefficients) then under appropriate conditions on the measurements, 0(S log N) are enough 
lfT9l . EOl . ||22l|. By developing similar tools, and applying them on the dynamics, rather than 
the signal directly, we show that with greatly-reduced control bandwidth, our algorithms perform 
close to the linear control bandwidth case. The performance is measured in terms of the distance 
between the recovered signal and the true signal. 

To the best of our knowledge, this is the first work in the space of wireless networks[] to 
exploit sparsity in the dynamics of the observed signal. As this is likely much more prevalent 
than sparsity in the actual trajectory of the signal (of course, if the trajectory is sparse, then so 
are the dynamics) we expect this high-level idea to find broad application. More concretely, the 

'Sparse changes in the dynamics of the signal have been used in other areas such as image processing to achieve video 
compression ITU . 
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main contributions in this paper are as follows: 

(1) A first (to the best of our knowledge) application of compressive sampling to reduce the 
control bandwidth in a cooperative sensing system. 

(2) A proof that path-loss matrices satisfy the null space property thereby allowing for efficient 
acquisition or sensing of the observed signal using standard convex programs such as i\- 
norm minimization and Lasso ll28ll . The proof technique is novel since path-loss matrices 
contain entries that have non-zero mean and are not independent, a scenario that has not 
been dealt with extensively in past research. 

(3) Simulation results that establish the competitive performance of our algorithm in compar- 
ison to the full control bandwidth case. 

The rest of this paper is organized as follows. In Section HH we introduce the system model for 
the cooperative sensing network under consideration. In Section HHl we discuss the compressive 
sampling algorithm that enables cooperative sensing using significantly-reduced control over- 
head. We establish the "goodness" of path-loss sensing matrices in Section [IV] The complete 
cooperative sensing algorithm is presented in Section |V] Simulation results establishing the 
competitive performance of the algorithm are contained in Section |VIJ Concluding remarks are 
made in Section IVIII 

Notation: Xij denotes element of matrix X while Xi denotes element i of vector x. (.) T is 
the transpose operator. For x E M. N , x^, A C {1,2,..., N} denotes the vector x restricted to 
the entries in A. For X e M mxn , X_4, AC {1,2,..., m} denotes the sub-matrix of X formed 
by the rows contained in A. Finally, 1 1 • | \ p is the p-norm operator on vectors. 

II. System model 

We consider a network with N s sensing nodes and a single fusion center operating in slotted- 
time. We introduce the necessary measurement and signal models below. 

Signal generation: We assume the signal of interest is being generated by multiple physical 
entities (e.g., TV transmitters) that are each dropped uniformly on a circle of radius r p cen- 
tered at the origin. There are a total of N p such signal generators that are located at points 
{(r p , Oj)}^! where 6j ~ £/[0,27r], Vj. The vector signal emitted at time t is denoted by 
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s(t) = [si(t) s 2 (t) . . . s Np (t)] T where s m (t) corresponds to physical entity or signal generator m 
and comes from some finite, discrete set. 

Spatial distribution of sensing nodes: There are N s , N s > N p sensing nodes that are placed 
on a collection of n c circles of radii {r S)1 , r S)2 , . . . , r s ric } where N s is such that ^* is even. Circle 
c contains ^ interfering receivers located at fixed points {(r S)C , m at are equally-spaced 

[Qi = ^pS i = 0, 1, . . . , ^ — 1 \ as shown in Fig. [IJ We note that this would roughly be the 
case when N s becomes large and the users are uniformly distributed on this collection of circles. 
The fusion center is located at some arbitrary point on the xy-plane. 

The above spatial distribution model affords us analytical tractability. In the simulations 
section, we show that the proposed algorithms work even under a more general spatial model 
where both the sensing and signal nodes are scattered uniformly on a regular square grid. For 
the sake of the analysis, we will also partition the sensing nodes according to the circle they 
belong to thus creating n c partitions {C\, C 2 , . . . , C nc } such that |J™! 1 d = {1, 2, . . . , N s } and 
Cj fl Cj = for i 7^ j. Within each circle, the nodes are numbered or ordered in diametrically 
opposite pairs as shown in FigfTJ, a labelling rule that is feasible since ^ is even. In other words, 
all pairs + 1) 6 Ci, j odd, will correspond to a pair of diametrically opposite nodes on 
circle Q. 




Fig. 1. Network with signal generators (not shown) uniformly distributed on the blue circle of radius r v . There are N s — 8 
sensing nodes in the network equally-divided across two circles (q = 2) of radii r s ,i and r S) 2 respectively. This gives rise to 
partitions Ci = {1, 2, 3, 4} and C2 = {5, 6, 7, 8}. The sensing nodes are equally-spaced on each circle as shown. 



Channel model: The distance between the n-th signal generator at location (r p , 6 n ) and m-th 
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sensing node at (r m ,8 m ) is given by 

d mn = V r m + r l- 2r m r n cos (9 m - 6 n ), 
which induces the following path-loss gain 

between the same. This model is an approximation of the free space path-loss model (with path- 
loss coefficient two) ETI . a choice that affords us analytical tractability while compromising 
very little on modeling accuracy. The composite channel gain between n-th signal generator at 
location (r n , 9 n ) and m-th sensing node at (r p , 9 m ) is given by 

h"mn(,t) 9mn (t) w^rrvn: (2) 

where g mn (t) ~ CN(0, 1) models small-scale Rayleigh fading. Distances do not change as a 
function of time since sensor nodes are deterministically placed and signal nodes are assumed 
fixed, once drawn from a uniform distribution. 

Sensor measurement model: We adopt the sensing model in ifTTI . lfT5l . where each sensor 
collects a set of samples over a window of size W time slots. In particular, at time t sensor 
node m receives the signal 

it) = ^2 h nm(t)s n (t) + V m (t), 
n 

where v m {t) is additive noise. Each sensor then forms the following measurement over a W- 
window of samples, 

(k+l)W 

z ^ k ) = W S Kit)? -n\v n {t)\\ (3) 

t=kW+l 

The sensors then send their measurements to the fusion center via orthogonal error-free control 
channels. The fusion center is charged with the task of recovering the observed signal {s n (t)}. 
In this work, we do not explicitly account for the time incurred in transmitting measurements to 
the fusion center and for subsequent actions such as data transmissions by the secondary network 
in the case of cognitive radio, as we are primarily interested in reducing the control bandwidth 
between the sensors and the fusion center. Thus, the time index t (and hence k) iterates only 
across sensing windows. 
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Sparse dynamics: Due to mismatches between the time-scale of the signal generators and the 
sensor sample collection period, we assume that the vector s(t) exhibits the following behavior. 
The signal s m (t) remains constant through the A;-th collection window kW + 1 < t < (k + 1)W 
and is denoted by s m (k). Furthermore, only a subset of signals change between the A;-th and the 
(k + l)-th collection window, i.e., |s(A; + 1) — s(A;)| is a sparse vector. This model is effective 
when the signal values change on a slower time-scale in relation to the sensing timeline. This is 
indeed what one typically expects for systems that exhibit a time-scale separation between the 
sensing network and the signal emitting process. Under the sparse dynamics assumption and for 
W sufficiently large, we invoke the Law of Large Numbers (LLN) and write the measurements 
in © as 

z{k)=As(k)+/3(k), (4) 

where A mn = E[|/i mn (t)| 2 ] from ©; j3 n (k) is a noise term that models the inaccuracies of 
the LLN over finite averaging windows. It is well-known from the Central Limit Theorem that 
E[|/3 n (&;)| 2 ] = O (jr) , Vn. The term (3(k) is referred to as LLN-noise through the remainder 
of this paper. Next, in the context of the measurement model in ©, we introduce the recovery 
algorithm that will be developed further through the course of this paper to exploit the sparse 
structure of the observed dynamics. 

Recovery algorithm: We assume that a subset Q C {1, 2, . . . , N s } of the measurements 
zq{k) are transmitted through \Q\ orthogonal control channels to the fusion center at the end 
of sensing window k. The fusion center must recover the signal s(k). The naive approach to 
recovering s(k) would be to transmit all the measured values, i.e., set Q = {1,2, .. . , N s } to 
the fusion centej 2 ). The above naive solution would consume a control bandwidth of 0(N S ). If 
the process s(k) is completely general, there is little that can be done to remedy this problem, 
and partial feedback (of only a subset of z) will necessarily result in degraded performance. 
However, as discussed in the introduction, for networks where only a subset of signal generators 
(e.g., TV transmitters) change their state between two adjacent sensing windows, we show that it 
is possible to reduce the control bandwidth. Using ideas from subset selection and compressive 
sampling, the next section considers how sparse dynamics can be exploited in order to achieve 
near-optimal performance while reducing the bandwidth consumption. 

2 The fusion center may recover the signal using standard tools such as least-squares. 
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III. Recovering dynamics through compressive sampling 

In this section, we propose a compressive sampling approach to efficiently recover the signal 
s(k). Before we discuss the recovery algorithm, we take a short diversion into the topic of 
compressive sampling. 

A. Compressive sampling 

The topic of compressive sampling has received tremendous interest in the recent years [[191 , 
[1201 . [|22l . The theory essentially states that one can recover sparse data exactly, given an under- 
determined system of equations. Specifically, the generic problem is the following: Given a 
signal x G M. N , one receives q « N linear, potentially noisy measurements: y = Mx + w. 
Here, M 6 ~R qxN encodes the measurement matrix and w G M. q denotes additive noise, usually 
of bounded norm. 

For a general vector x G M. N , N independent measurements are required to hope to reconstruct 
x. When q < N, the problem therefore is under-determined. If x is sparse, however, in some 
settings the problem is no longer under-determined, and can be solved exactly by considering 
standard convex programs such as 

arg min : I |Mx — y| I, — £||x||i for some chosen f > 0, 

X6RP " 112 " " 

in the noisy case and 

arg min : I |x| N 

(5) 

Mx = y. 



in the noiseless case. The former is the so-called Lasso 11281 formulation for model selection 
(subset selection) while the latter 1'i-norm minimization problem is called Basis Pursuit. Many 
such results have appeared in the literature, e.g., |fT9l , 11201 . [|22l . 0711 under the umbrella of 
compressive sampling^. Indeed, the results are attractive from an algorithmic perspective since 
the convex relaxation is easily solvable, with computation time that scales gracefully as the 
size of the problem increases, allowing the efficient solution of very large problems. The above 
convex programs succeed as long as the linear equations, or measurements, satisfy a property 
called Null Space Property (NSP), which essentially amount to the statement that there are no 



3 We do not use the terminology "compressed sensing" to avoid confusion with "cooperative sensing". 



10 



very sparse vectors in the null- space of the measurement matrix M. The theoretical connections 
between Lasso and Basis Pursuit have been well- analyzed by authors such as Tropp 11381 . 

B. The recovery algorithm 

Returning to our problem, we define As(fc) = s(k) — s(k — 1) and apply the model selection 
paradigm outlined above, to the dynamics vector As(k) rather than the signal vector itself. We 
can assume that at some initial time k , s(k ) is known. At time k, we "query" a subset of sensors 
Q and receive measurements z Q (k) = A Q s(k) + (3(k). We can then construct the difference in 
measurements 

Az Q (k) = z Q (k) - z Q (k - 1) 

= A e [s(A;)-s(A;-l)]+/3(fc)-/3(fc-l) 
= A e As(fc)+/3(fc)-/3(fc-l). 

Since the left hand side, Azg(fc), is known, this falls precisely into the sparse recovery paradigm 
developed above. More concretely, let Q be the subset of queried users with size \Q\ = q, q even. 
Then, Q is chosen according to the following algorithm: 

Note that the output query set is not unique since the choice of node-pairs is left open. Let Q 

Algorithm 1 Protocol to choose query set Q of size \Q\ — q 
i: Set Q = and i = 0. 

2: while % + 1 < | do 

3: Choose any pair of diametrically opposite receivers + 1) from circle % + 1, i.e., 
j e C i+ i, j odd. 

3: Q = QU{j,J + 1}. 

3: SetC i+1 =C l+1 \{j,J + 1} 

4: Increment i = (i + l)(mod n c ). 

5-. end while 



be the set of all possible output query sets from Algorithm Q] above. For the example in Fig. [T] 
one possible output query set for q = 4 is Q = {1, 2, 7, 8}. It is necessary for a query set Q to be 
selected in this way for the sake of analytical tractability. As with the spatial distribution model, 
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we adopt a more general querying model in our simulations section. Under such a querying 
model, it is of immediate interest to determine the smallest query size q (or control bandwidth) 
that the system requires in order to recover s(k) reliably using 

argmin xeR iv P : ||A Q x - Az Q (fc)||^ - f ||x||i. (6) 

As mentioned earlier, in this work, we do not consider the number of bits required to commu- 
nicate zq(/c) reliably to the fusion center as we are interested primarily in the scaling behavior 
of the control bandwidth. 

Compressive sampling theory states that recovery of any S'-sparse vector through Lasso or 
Basis Pursuit is possible in a noiseless setting if and only if the sensing matrix Aq, Q E Q 
satisfies the NSP ll32ll of order S. This property is defined in the next section. We note that in our 
application, the sensing matrix is provided by the channel as opposed to traditional compressive 
sampling where the designer is allowed to choose a convenient sensing mechanism. In the next 
section, we present the main result of this paper, which states path-loss matrices Aq, Q G Q 
make for good sensing matrices in the noiseless case, i.e., when W — > oo. Recovery results that 
are proved for Basis Pursuit, which, in the absence of noise, carry over to Lasso 11381 . While 
the theory is developed for the noiseless case with W — > oo, our simulations consider finite 
averaging windows and demonstrate that sparse recovery still remains effective in this setting. 

IV. NSP OF PATH-LOSS MATRICES 

In this section, we establish that path-loss matrices Aq, Q G Q satisfy the Null Space Property 
(which will be defined shortly) when the control bandwidth obeys q = 0(S\og N p ). Lemma HI 
Lemma [5] along with Theorem @] constitute the main results in this section. 

A. Preliminaries 

Let the support set of x be denoted by S. A vector x is S'-sparse if |«S| < S. We define the 
null space property from Gribonval et al. Il34ll . Given a matrix M, let A/"(M) denote its null 
space. 

Definition (Null space Property): A matrix M satisfies the null space property of order S if for 
all subsets 5C{1,2,...,JV} with |«S| < S, the following holds 

||v 5 ||i < ||v 5 c||i, Vv e Af(M) \ 0. 
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where S c = {1,2,..., N} \ S. Based on this property, the following recovery result [34] has 
appeared both implicitly and explicitly in works such as 11321 . 11351 . 

Theorem 1: [34] Let M £ W xN . Then, any ^-sparse vector may be recovered by solving 
© iff M satisfies the NSP of order S. 

□ 

The NSP is typically quite difficult to prove directly leading to the development of sufficient 
conditions that are easier to establish. One such sufficient condition is the restricted isometry 
property []37l that has become quite popular in recent years and is defined below. 
Definition (Restricted Isometry Property): A q x N matrix M satisfies the Restricted Isometry 
Property (RIP) of order p if there exists e p (M) £ (0, 1) such that for all x £ M. N , 

(1 - e p (M))||x r ||2 < ||x^(M T ) r |l2 < (1 + e p (M))||x T ||2 (7) 

holds for all sets T with |T| <p. Here, (M. T )j- is the sub-matrix of M T formed by rows in T. 

Here, e p (M) is called the restricted isometric constant of M. The RIP essentially requires that 
all q x |T| sub-matrices of M be well-conditioned. Under such a conditioning, perfect recovery 
of x is possible as stated in the following theorem. 

Theorem 2: flH, flU Let M £ WL qxN . If M satisfies the RIP with e 25 (M) < « 
0.4531, then every 5-sparse vector x £ M. N is the solution to the £i-norm minimization problem 
in ©. 

□ 



Thus, the RIP with a sufficiently small constant immediately implies the NSP in the context of 
£i-recovery. The approach we use to prove "goodness" of path-loss matrices Aq is motivated by 
the following observation. In general, the null space of a product of two matrices NM contains 
the null space of M and therefore if NM satisfies the NSP, so does M. This allows us to study 
the class of linearly-processed path-loss matrices Aq = B WAq where 



W = diag{JJ .^J}, J 



f times 



(8) 
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and diag{-} is the standard block-diagonal operator; 




support {±1}. We focus our attention on establishing the recovery properties of Aq rather than 
Aq. We show that Aq satisfies the RIP with q = 0(Slog N p ) measurements and hence the NSP. 
The transformation W essentially subtracts rows of Aq corresponding to diametrically opposite 
pairs of sensors in the same partition. Thus, the dimension of G is still q x q. The transformation 
B weights and adds adjacent rows of G. 

According to our spatial distribution model, since the positions of the sensing nodes are 
fixed, the columns of Aq become stochastically independent since each interfering transmitter 
is independently thrown. At this point, we will rely heavily on recent results from Vershyin 
[|39l and Adamcyzk et al. [|40l that deal with sensing matrices containing independent columns. 
Before we reproduce the RIP result 11391 . ||40| for matrices with independent columns, we present 
a primer on sub-gaussian and sub-exponential random variables along with some useful results 
from non-asymptotic matrix theory. 

B. Useful concentration inequalities 

We refer the reader to the tutorial paper by Vershynin ll39l for a great introduction to non- 
asymptotic matrix theory. Lemmas |T]{3] below are well-known past results that are summarized 
in this paper ||3~9~1 . The proofs are not reproduced due to lack of space. 



Lemma 1: Let u be random variable. The following properties are equivalent with parameters 
Ki > 0, i — 1, 2, 3, 4, differing from each other by at most an absolute constant factor. 

(i) Tails: Pr(|u| > 7) < exp(l - for all 7 > 0, 

(ii) Moments: (E [\u\ p ])p < K 2y /p for all p > 1, 



Moreover, if K[u] = then properties (i)-(iii) are also equivalent to the following property: 
(iv) Moment generating function: E [exp (try)] < exp(7 2 i^ 4 ) for all 7 G M. 



A random variable that satisfies the above property is called a sub-gaussian random variable. 




□ 
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Such random variables are often characterized by the ^-normj, which is defined as 

(E\\u\ p ])p 
IHk = sup p > 1 — — . 

It follows that if the ^ 2 -norm of u is finite, then u is a sub-gaussian random variable with 
| \u\ = K 2 . This is in fact the case for bounded random variables with symmetric distributions. 



Lemma 2: Let u be a symmetrically distributed, bounded random variable with \u\ < M, M > 
0. Then, u is a sub-gaussian random variable with ||m||^ 2 < cM 2 , c > 0. 

□ 

In higher dimensions, a random vector u of dimension N is called sub-gaussian if u T x is 
sub-gaussian for every x G M. N . 

Lemma 3: Let {«i}f =1 be a collection of independent, zero-mean, sub-gaussian random vari- 
ables. Then, u is a sub-gaussian random vector with ||u||^, 2 = C maxj ||wi||^ 2 for some C > 0. 

□ 

We are now ready to prove the RIP (hence NSP) for matrix Aq. Before we move on to this 
task, we require one more definition. A random vector m of dimension N is called isotropic if 
E[|m T x| 2 ] = ||x|| 2 for all x G R N . 

C. NSP of linearly-processed path-loss matrices Aq 

We reproduce the recent RIP result ||39ll , HU1 concerning matrices with independent columns. 
We refer the reader to |[39ll . Il40l for the proof. 

Theorem 3: Let M = [mi m 2 . . . m N ] be a q x N random matrix whose columns are 
independent, isotropic and sub-gaussian with ip m ax,m = maxj 6 /i ,...,jv} H m ilU 2 - Furthermore, 
let the columns satisfy 1 1 1 1 2 = q almost surely for all i. Then, the normalized matrix is 
such that if q > C^ max m e~ 2 S\og j f or some £ > 0, then e p (^ M ) < £ witn probability 

at least 1 -2exp(-c^ max m e 2 q). Here, c^ max m and C^ max m are positive scalars that depend only 
the worst-case sub-gaussian norm ip max ,m- 

□ 

As mentioned earlier, the channel matrix A = [Ai A 2 . . . A^p] contains independent columns 



4 Alternate definitions of this norm have been adopted (such as in 1401 ) that are all equivalent to within a constant factor. 
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since the positions of the sensing nodes are fixed. However, each column contains entries that 
are not centered, not isotropic and that are highly coupled. This is because all entries in Aj 
are now completely determined by the position of signal generator i. For this reason, it is not 
immediately clear whether the columns are sub-gaussian. 

To finally prove the NSP of Aq, Q E Q, we show that after left-multiplication by matrices 
B and G, the matrix Aq satisfies the sufficient conditions in Theorem [3] The following lemmas 
and theorem constitute the main theoretical results of this paper. 

Lemma 4: The matrix Aq = BWAq, Q E Q of size q x N contains independent, isotropic, 
centered, sub-gaussian columns. 

Proof: See Appendix lAl □ 

Lemma 5: For matrix Aq = BWAq, Q E Q of size q x N, we have that ||ai|| 2 = q almost 
surely. 

Proof: See Appendix [Bj □ 

Theorem 4: Aq, Q E Q satisfies NSP of order S almost surely when q = 0(SlogN p ). 
Proof: The result follows from Lemma H Lemma [5] and Theorem [3] □ 



We discuss in the next section, how this sparse recovery algorithm is incorporated into the 
development of a complete cooperative sensing algorithm. In particular, since we are estimating 
dynamics rather than the signal itself and due to the presence of LLN-noise, the real possibility 
of error propagation arises. This has not been heretofore addressed in the literature, to the best 
of our knowledge. In the next section, the algorithm introduces explicit steps to control this. The 
simulations in Section [VI] demonstrate the effectiveness of our approach. 

V. Complete cooperative sensing algorithm 

The essential conclusion of the previous section is that the signal vector s(k) can be recovered 
by acquiring 0(S\ogN p ) noise-free measurements if the dynamics of s(k) are sparse. 

The complete sensing algorithm is presented in Algorithm [2] Recall that we do not consider 
the actions taken post-sensing and the time taken (equiv. bandwidth per control channel) to 
transmit the measurements to the fusion center as we are only interested in scaling properties 
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in this work. The algorithm operates in three modes; in dynamics mode, the algorithm retrieves 
signals dynamics incurring logarithmic control bandwidth and then estimates the recovery error. 
If the estimated recovery error is too high, the algorithm switches to partial reset mode in the 
next sensing window or slot where the fusion center retrieves the signal directly while incurring 
linear bandwidth. Following partial reset, the system returns to dynamics mode. This operation 
is overlaid by a periodic (with a pre-set period) reset mode where again the fusion center again 
retrieves the signal directly incurring linear bandwidth. 

In a typical sensing slot, the Lasso procedure in Step [8] recovers As(fc). The recovery is almost 
surely exact when As(fc) is sparse and when W is sufficiently large. However, as this might 
not always be the case in reality, error propagation is a critical aspect that needs to be managed 
when operating on the dynamics of the signal. The following remarks describe the steps taken 
by the algorithm to address the same: 

• The system selects a parameter T err that represents the number of sensing slots after which 
the fusion center must query all sensors and recover the absolute signal as opposed to the 
dynamics. 

• During typical operation, i.e., when /c(mod T err ) ^ 0, we estimate the accuracy of s(k) in 
Step [TOl immediately following compressed sampling in Step [9] The technique is commonly 
referred to as cross-validation ||4TI . If the estimated error is larger than the threshold e, we 
query all remaining N s — (q + q') sensors and recover the absolute signal in the next time 
slot. 

• Once more, due to the sources of noise described above, the vector u(k) retrieved by Lasso 
in (flOl) and ([12]) will almost never be truly sparse. Thus, we introduce a thresholding function 
to recover a discrete signal. 

VI. Numerical results 

In this section, we evaluate the performance of the proposed sensing algorithm using numerical 
experiments. 

A. Simulation setup 

The simulation setting is described in the following. 
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Algorithm 2 Sensing algorithm 



1: Primary inputs: sensing window sizes W, number of sensors N s , number of signal emitters N p , sparsity factor 

S, error reset period T err > 
2: Initialize mode[0] = "dynamics" 
3: for each k = 1, 2, . . . , do 
4: if fc(mod T err ) then 

5: During time slots kW < t + 1 < (k + 1)W, each sensor samples the signal of interest W times and 

averages the same to obtain z n (k) as per @. 
6: if mode[k-l] = "dynamics" then 

7: The fusion center queries q = Q(S\ogN p ) sensors indexed by set Q according to Algorithm [T] and 

computes the change in received signal levels {Azg(fc)}. 
8: Dynamics mode: Once acquired, the fusion center estimates the change in measurements As(k) = 

s(k) — s(k — 1) by solving the following optimization problem: 

u(fc)=arg min : ||A Q x - Az Q (fc)|| - fllx^, (10) 

and setting 

Ai(fc) = [#(!!(*)) ^(U(fc)) ^ p ( U (fc))] T , 

where <\>\ : M. p — > R is a suitably-chosen thresholding function with v that is used to recover discrete 
sparse signals in noisy environments. 
9: The fusion center sets s{k) = s(k — 1) + As(k). 

10: The fusion center then queries and acquires extra measurements from q' additional sensors Q' C 

{1, 2, . . . , N s } \ Q, \Q'\ = q'. Then, the fusion center computes 

e=^\\A Q ,~s(k)-z Q ,(k)\\ 1 , (11) 



which essentially represents the average recovery error on the extra measurements, 
if e < e, i.e., the estimate is accurate, then 

Set s(fc) to be the final estimate of the signal. 

Set mode[k] = "dynamics" 
else 

Set mode[k] = "partial reset" 
end if 

else if mode[k-l] = "partial reset" then 

Partial Reset Mode: The fusion center queries all N s sensors and uses the measurements to directly 
recover the final estimate s(k) by solving 

u(fc)=argmin : ||A Q x — za( fe )lll — ^I! x l|i ( 12 ) 

and again, setting 

s(fc) = k(u(fc)) #(u(fc)) p N (u(k)) 1 



19: Set mode[k] = "dynamics" 

20: end if 

21: else if fc(mod T err ) = then 

22: During time slots kW < t < (k + 1)W, each sensor forms z n {k) as per (0]). 

23: Full Reset Mode: The fusion center then queries all sensors and uses the N s measurements to directly 

recover the final estimate s(k) by solving ( fT2b . 
24: Set mode[k] = "dynamics". 

25: end if 
26: end for 
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Network geography - We consider a square region with side length 1 kilometer. We partition 
this square area into N = 49 micro-cells, each micro-cell is 143m x 143m square area, which 
is sufficiently small to ensure at most one signal generator or sensor per micro-cell. We index 
the micro-cells by i (1 < % < N), counted column-wise as shown in Figure [21 Each micro-cell 
contains one potential signal generator that is located at the center of the square micro-cell. The 
sensors are positioned one per micro-cell, i.e, N s = N = 49. 
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Fig. 2. Illustrative example of grid model used in simulations with N — 25 and 7V S = 4. 

Signal dimension and channel state knowledge - By virtue of the above geographical model 
where we have "discretized" the locations of the sensors and signal generators, the fusion center 
may now pre-compute A based on the grid. This removes the need for real-time knowledge about 
the positions of the signal generators. However, this comes at a price, namely, an increase in the 
dimension of the signal from N p to N. Thus, the value s m (k), m — 1, 2, . . . , N, is non-zero if 
a signal generator is present and active at micro-cell m on the grid. 

Signal activity - In our simulations, we assume that the signals are binary-valued, i.e., s m (k) e 
{0, P}, P > 0, Vm, k. The activity of each signal generator or micro-cell s m (k) is modeled as 
a Markovian ON-OFF process where the transition probability from to P is 770 = rj = 0.1 
and from P to is r]i = 77 = 0.1. The probability that the signal changes state from sensing 
window to the next can be computed as = 2 (o 1+0 i') 1 = ^-l- Thus, on average 10% or 

roughly S ave = 49 x 0.1 ~ 5 signals change their state, which reflects the sparse dynamics. 
As per our system model, we assume that P is known and thus, we are only interested in 
estimating the support at time k. Henceforth, the true support at time k be given by V(k) = 
{m : m = 1,2, . . . , N, s m (k) = P}, while the estimated support by Algorithm [2] is given by 
V(k) = {m : m = 1, 2, . . . , N, s m {k) = P}. 
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Post-recovery thresholding - We use a simple per-element thresholding function 



if 



signfxjP, \xA > v 
» = { y ! ' ' - (13) 
0, else 

to recover the actual sparse entries themselves. Note that we have to account for negative values 
for Xi since u(k) represents differences in ( fTOl ). 

Performance metric - To measure the signal recovery accuracy of our algorithm, we use a 
set distance. Given two sets A,BC {1,2,..., N}, the measure of distance that we use is 

d(A,B)=l- — [\AnB\ + \A c nB c \}. (14) 

In the context of our sensing problem, i.e., setting A = V(k) and B = V(k), we see that d(-) 
essentially counts the fraction of elements that are incorrectly recovered at time k. Now since 
we are dealing with a stochastic, time-varying system, we are interested in the long-term average 
distance that is given by 

1 T 

-J2d(V(k),V(k)), (15) 
fc=i 

where T is the length of our simulation counted as the total number of sensing windows. 

Other algorithm parameters - The number of measurements for recovery and error mea- 
surement are set to q = S ave logiV = log N = 77 log N = 29 and q' = 5 with associated 
query sets Q = {1, 2, . . . , k} and Q' = {k + 1, k + q'}. The regularization parameter is set to 
£ = 0.01. The remaining parameters such as W,r),e and T err will be varied in our experiments. 

Control bandwidth - Note that the parameters q, q', e, T err and N have a direct bearing on 
the bandwidth consumed by our algorithm. To analyze the bandwidth consumed in Algorithm [2l 
we basically count the fraction of time algorithm spends in each of the modes in Steps [8l [T8] and 
l23l For a total communication duration of T sensing windows or slots, i.e., k — 1, 2, . . . , T, let 
the time spent in dynamics, semi-reset and reset mode be given by Td, T s __ r and T r respectively. 
Then, the bandwidth consumed, is given by 



B = ZL( q + q >) + N< T — + T " 



ZL( q + q >) + N(l-ZL 



(16) 



and the savings over full (linear) control bandwidth can be expressed as 1 — J|. The control 
bandwidth, in addition to the distance measure in (fT4l) . is an important evaluation metric, one 
that forms the core motivation of this paper. 
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Through the remainder of this section, we explore the many interesting trade-offs that are 
inherent in our proposed sensing design in the context of the above parameters. 

B. Simulation results 

Having described our simulation setup in detail, we now proceed to analyze the performance 
of our algorithm. Since the baseline naive algorithm in itself does not achieve perfect recovery 
due to LLN noise, we adopt the following methodology for comparison. We fix a target average 
distance and determine the minimum sensing window size W m i n needed to achieve this target 
when consuming full linear bandwidth, i.e., with q = N. Note that as the sensing window size 
grows, the average distance can be made arbitrarily small thereby approaching perfect recovery. 
We show that with careful tuning, Algorithm [2] consumes considerably lesser bandwidth in 
achieving the same target average distance while maintaining the same window size of W min . 

First, in order to determine W min , we analyze the average distance of the naive algorithm 
versus W, i.e, set T err = 1 in Algorithm [2l which means that we always retrieves absolutes. We 
fix the target average distance to be 0.1. In Fig. [3l we plot the average distance of the naive 
algorithm over the range W E {50, . . . , 1000}. From Fig. [3l we see that for a target average 
distance of 0.1, we need 2000 samples to form (HJ). We therefore set W m i n = 2000 as the baseline 
window size for the naive algorithm. 




Fig. 3. Average distance of the full linear bandwidth algorithm versus W, i.e, T err — 1. 
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Having identified W min , next we proceed to tune Algorithm [2] with the goal of outperforming 
the naive algorithm. In particular, we study the accuracy-bandwidth performance of our approach 
while varying the error threshold e, the reset period T err and the support threshold v over the 
space (e,T err ,u) E {0.1P, 0.5P, P, 2P, 3P} x {P, 2P, 3P, 4P, 5P, 8P} x {0.25P, 0.5P, 0.75P}. 
Note that this amounts a total of 5 x 6 x 3 = 90 configurations. The window size is set to 
W = W m in = 2000 for a fair comparison with the baseline case. The number of measurements 
in dynamics mode is set to q = \S ave logN~\ = 29. The number of extra measurements for 
error control is set to q' = 5 ~ 0.15g. The average distance and bandwidth saved under each 
configuration are plotted in Figs |4(a)| and Figs |4(b)| respectively. The x-axis here is a configuration 
index, which essentially maps to a (e/P, T err /P, v/P) -triplet as specified in Table [fl From 
Table U we see that configurations 71 to 90 meet the average distance target of 0.1 and in 
fact, configurations 71 to 80 do so while saving up to 15% in control bandwidth. 




10 20 30 40 50 60 70 
Configuration index 




10 20 30 40 50 60 70 80 90 
Configuration index 



(a) (b) 
Fig. 4. (a) Average distance versus configurations, (b) Bandwidth savings versus configurations. 



Thus, by exploiting the structure in the emitted signal, we are able to achieve significant savings 
in control bandwidth over the naive algorithm. Recognizing that the number of parameters that 
control the algorithm is quite large, we believe that with more extensive tuning/optimization of 
the various components, e.g., error measurement budget q', further savings might be possible, 
especially as the scale or dimension of the problem grows. Finally, we note that the proposed 
compressive sampling approach has been applied to (even) more general spatial distribution 
models such as one where both the sensor and signal nodes are thrown uniformly at random on 
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a square area, in the absence of a grid. The results are as promising but not presented here due 
to lack of space. 

VII. Conclusions 

In this paper, we apply ideas from the field of compressive sensing to reduce the control 
bandwidth in cooperative sensing systems. We are particularly interested in settings where 
the dynamics of the signals being observed are sparse. This is indeed what one typically 
expects for systems that exhibit a time-scale separation between the sensing network and the 
signal emitting process. We prove theoretical RIP-/NSP-based recovery results for path-loss 
matrices in networks with circular geometries. We then develop a complete cooperative sensing 
algorithm that addresses key issues such as error propagation. As the proposed algorithm requires 
many input parameters, we extensively tune these parameters in our numerical experiments. We 
demonstrate that our approach can provide significant bandwidth savings over the naive linear 
bandwidth case under more general spatial distributions. Extending the theory to cover these 
general spatial distributions remains an open problem that deserves attention. 
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Recall that A c . represent the sub-matrix of Aq containing the rows specified in Cj. Then, the 
entries {Ay } of sub-matrix A Cj are identically distributed. This follows from the observation that 
any two sensor nodes on the circle of radius r s>i will perceive the same distribution of signal 
nodes since the latter nodes are distributed on a circle. Recall that the transformation G = WAq 
essentially subtracts rows of Aq corresponding to diametrically opposite users on each circle. 

Thus, the columns of G are independent and all entries are centered. The next step is to show 
that the entries = Ay — Xa+i)j are symmetric. This is does not follow immediately from 
the fact that and A(j + i)j are identically distributed when indices i and i + 1 come from the 
same partition since they are not independent. Hence, we appeal to the concept of exchangeable 
random variables, which is defined below for the special case of a pair of random variables. 

Definition (Exchangeability): Two random variables X and Y are called exchangeable if their 
joint cumulative distribution function (cdf) is symmetric, i.e. if F x ,Y(x,y) — Fx,y(y,x). 

It is known the difference of two identically distributed, exchangeable random variables is 
indeed symmetric ll43l . Thus, we only need to establish this fact for the case when A^ and Xu+i)j 
come from the same circle, say C c , since the definition of G in ([8]) precludes any other possibility. 
To establish that A^- and are exchangeable, we compute the joint cdf F Xlj (x, y), 

when i and i + 1 come from the same circle c, as follows 



Appendix A 



Proof of Lemma [4] 




0, y) = Pr(A y < x, < y) 

= Pr((% >\-K, d(i +X )j > 



— y 



i 



K) 




2jt 



where we have introduced the function 



d 2 (8 m , 9 n ) 



r li + r l ~ 2r p r s>1 cos (8 m -9 n ), i = 1, . . . , a(k), 



(17) 
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for notational convenience through the remainder of the proof. The set 71(9 j, 9 iy x, y) is defined as 
7Z(9„9 t ,x,y) = Uj :d(9„9 t ) >[l-K] + ,d 2 > [j-^] + J and V ( A ) denotes 

the volume of the set A. The volume of set 71(0 j, 9 U x, y) can be expressed as a sum of the 
volumes corresponding to smaller sets. To that effect, we define 

7l s (0,a) = jfy : d(0j,0) < ^-K J 
and can thus write V(7l(0 j ,0 i , x, y)) = 2nr p - [v(7Z s (9 h x)) + V (k s + |f , yjj - V (7Z S (0 U x) 

From this characterization, we can immediately conclude that V(7Z(r p ,9j,r sc ,9i,x,y)) = 
V(7Z(9j, 9 iy y, x)), which gives our the desired result that A^ and \u+i)j are exchangeable 
implying that the entries g^ = A^ — are symmetric for all (i,j). It is clear from the 

above arguments that exchangeabilty essentially follows due to the uniform distribution of the 
signal emitters. 

The symmetry of gij is crucial for our next step where we argue that columns of Aq = BG 

remain independent. By definition, = , 1 --^%giy Thus, to prove that the columns of 

y Var{gii} 

Aq = [ai a.2 . . . ajVp] are independent, we need to show that _L for any arbitrary i,j 
and k ^ j. Since the Bernoulli random variables are independent across rows and since g^ _L 
g m k, k 7^ j, m ^ i, we clearly have that aij _L a m k, k ^ j, m ^ i. Thus, we only need to 
establish that _L for k ^ j or equivalently that fiigij _L fiigih- But this follows from 
the symmetry of g^ which means that knowledge of ^%gij reveals no information about the 
random variable Thus, the columns of Aq are independent. In addition, there are identically 
distributed and hence we can now focus on studying the properties of column a 1 without loss 
of generality. 

The random vector ai is isotropic since Efaf-J = j^^E [g^] = 1 and 

E[a a a kl ] = —2=—2==E [f3, t f3 k ] E \gn9ki] = 

for % ^ k, implying that E[|afx| 2 ] = Ei E [ a ?iR 2 + Y,i^k E i a n a ki] x i x k = ||x|| 2 . 

To prove that a x is sub-gaussian, we first condition on the position of the first signal emitter 
(r p , 9\). This will allow us to apply Lemma [Hand Lemma [2] since ga is now completely known 
thereby making a a a collection of independent random variables. The elements an are symmetric 
and bounded with |a 7l | < 1=== when conditioned on position (r„,6-i). 
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Lemma 6: There exists M p > 0, p = 0, 2, . . . , q — 1 such that 

E[<&] = M p , for < i < (p + 1)^, Vfc. (18) 

Proof: By definition, the distribution of gn for < i < (p + 1)— depends on the distance 
between the corresponding diametrically opposite sensors on circle (p + 1) along with the 
distribution of the first interfering user. Since this distance always remains the same independent 
of % and k, the result follows. □ 

Hence, by Lemma [2] and since |<7a| < 1, an is sub-gaussian with | \an\ \^ 2 \{r P ,6i) < ^p. where 
— min p= i ^ q Mp. Here, we have introduced notation 1 1 • | \( rp ,6i) to indicate explicitly that we 
have conditioned on the location of the first interfering user. Now, from Lemma [U we conclude 
that ai is a sub-gaussian vector when conditioned on the location of the first signal node with 
| |ai ||,/, 2 \(r p ,6i) < 2C, C > 0. However, since the sub-gaussian norm ||ai| |^ 2 |(r p ,fli) com P ute d 
above is independent of (r p , 8i), this implies that ai is a sub-gaussian vector with ||ai||^, 2 < 2C. 
This can be seen by applying the Law of Total Probability to the definition of sub-gaussianity 
in Lemma [Q □ 

Appendix B 
Proof of Lemma [5] 

To show almost-sure convergence of the norm of sl 1 , we need to prove that 3c* > such 
that Pr (lim^oo | \ £? = i 4 " c * | = °) = Pr ( lim fe^oo | \ Eti & ~ c* | = o) = 1 where the 
probability is computed over the random location (r p ,9), 9 ~ Z7[0, 27r]. We instead prove the 
following more general statement that 3c* > such that 



lim 



k— >oo 



9tx 



for all (r p ,9), 9 E [0,2tt]. (19) 



Note that (fT9l is a completely deterministic convergence statement in contrast to the earlier 
probabilistic statement. From the proposed query protocol in Algorithm [H we see that the number 
of sensors selected for feedback is a monotonically increasing function of k for all circles. 
This mean that we can study the convergence of the norm for any one circle (a sub-vector 
of ai) and draw conclusions about the norm concentration of the entire vector ai. Consider 
the first partition of sensors on circle C\ and let the size of this partition be a(k). Then, we 
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see that a(k) — > oo as k — >■ oo and hence, we can shift our focus to proving the property 



lim 



k— >oo 



a(k) 2^,1=1 



c* = 0, V(r p ,0i), 6»i e [0,2n]. In fact, we show that c* = 1, i.e, 



limt. 



k— >oo 



— V 



a(fc) 2 



- 1 



0, V(r p ,0i), 0i G [0,2tt] 



(20) 



Recall that the squared-distance between a signal generator located at (r p ,0i) and sensor node 
i on C\ is given by d 2 (6i,9i) where d(-) is defined in (fT71)and where 0, = ^rp. Now let 



f(x,y) 



x, y > where K is given in the definition (Q~|). Then, can 



K+x K+y 

equivalently be expressed as g\ = f (d 2 (0j,0i) ,d 2 (0; + §,0i))- Now, we proceed to prove 
(|2Q|) by observing that 

1 V^a(fe) 2 



lim fe _ 



(21) 



= lim^oosfe Efi? / (*» *i) , d 2 + ^Mi)) 

= h C f ( d2 *i) ,<* 2 (0* + f , *i)) ^- 
This claim follows from that fact that the expression on the left is essentially the Riemann sum 
of the integral on the right. This means that for any given e > 0, we can find k £ such that when 



k > k e , we have that 
calculate 



a (fc e ) Z-fi=l i/il 2vr JO ^ I" 



1 r 2?r 

2~ 



7T 



<f 0, + -,0 d0. 



< e. Next, we 



(22) 



By substituting (fT7l) in (|2TT) and (T22l) . we see that the variable of integration may be switched 
and hence 

= f{d 2 {6 u 6),d 2 {6 l + ^6))d6 



d 2 (e i + i,e))dBi 



(23) 



which implies that 



a (fc e ) Z^i=i »ti 2tt Jo J I" 



jd 2 (fl. + jE^))^. 



(24) 



We can then establish convergence through 



1_ v^a(,fe 
ffel Z-a=l 



«(*) 9? 



«(*) 



1 



_J_ S^a(k) 2 _ E r„2 1 



*(*) 



< 



e4] £ when / ' :>/ ' • 



The result follows since E^] is bounded below by Lemma [6] 



Config. 


[e/P, T err /P, WP)-triplet 


Distance 


Bandwidth 


index 






savings 


1 


(3,0.25,0.8) 


0.31 


0.2634 


2 


(2,0.25,0.8) 


0.2838 


0.2499 


3 


(3,0.25,0.5) 


0.2588 


0.2432 


4 


(3,0.75,0.8) 


0.253 


0.2597 


5 


(3,0.5,0.8) 


0.2454 


0.2667 


6 


(2,0.25,0.5) 


0.2442 


0.2338 


7 


(2,0.5,0.8) 


0.2336 


0.2584 


8 


(3,0.25,0.4) 


0.2327 


0.2285 


9 


(2,0.25,0.4) 


0.2244 


0.221 


10 


(2,0.75,0.8) 


0.2227 


0.2381 


11 


(1,0.25,0.8) 


0.2091 


0.2025 


12 


(3,0.75,0.5) 


0.2076 


0.2383 


13 


(3,0.25,0.3) 


0.1999 


0.2037 


14 


(1,0.25,0.5) 


0.1941 


0.1904 


15 


(2,0.25,0.3) 


0.1933 


0.1997 


16 


(3,0.75,0.4) 


0.1909 


0.2233 


17 


(2,0.75,0.5) 


0.1898 


0.218 


18 


(3,0.5,0.5) 


0.1888 


0.2445 


19 


(1,0.25,0.4) 


0.1875 


0.1836 


20 


(2,0.5,0.5) 


0.1848 


0.2405 


21 


(2,0.75,0.4) 


0.175 


0.2064 


22 


(1,0.75,0.8) 


0.1724 


0.1863 


23 


(3,0.75,0.3) 


0.1693 


0.1988 


24 


(1,0.25,0.3) 


0.1683 


0.1663 


25 


(1,0.5,0.8) 


0.1683 


0.2207 


26 


(3,0.5,0.4) 


0.1635 


0.2295 


27 


(2,0.75,0.3) 


0.1628 


0.1802 


28 


(2,0.5,0.4) 


0.1603 


0.2261 


29 


(1,0.75,0.5) 


0.1597 


0.1694 


30 


(1,0.75,0.4) 


0.1568 


0.1733 


31 


(0.5,0.25,0.4) 


0.1562 


0.1552 


32 


(0.5,0.25,0.8) 


0.1562 


0.1579 


33 


(0.1,0.25,0.4) 


0.1514 


0.1531 


34 


(1,0.25,0.2) 


0.1511 


0.1531 


35 


(0.1,0.25,0.2) 


0.151 


0.1531 


36 


(3,0.25,0.2) 


0.151 


0.1531 


37 


(0.1,0.75,0.8) 


0.1508 


0.1531 


38 


(0.1,0.25,0.8) 


0.1506 


0.1531 


39 


(2,0.25,0.2) 


0.1497 


0.1531 


40 


(0.5,0.25,0.2) 


0.1494 


0.1531 


41 


(0.5,0.75,0.8) 


0.1491 


0.1592 


42 


(0.1,0.75,0.4) 


0.1474 


0.1531 


43 


(0.5,0.75,0.2) 


0.1472 


0.1531 


44 


(1,0.5,0.5) 


0.1463 


0.2063 


45 


(2,0.75,0.2) 


0.1461 


0.1531 


46 


(1,0.75,0.2) 


0.146 


0.1531 


47 


(0.5,0.25,0.5) 


0.1453 


0.1393 


48 


(0.1,0.75,0.2) 


0.1448 


0.1531 


49 


(3,0.75,0.2) 


0.1448 


0.1531 


50 


(0.5,0.75,0.4) 


0.1446 


0.1569 


51 


(1,0.75,0.3) 


0.1444 


0.1422 


52 


(0.5,0.75,0.5) 


0.1434 


0.1374 


53 


(0.5,0.75,0.3) 


0.1392 


0.1151 


54 


(0.1,0.75,0.5) 


0.1391 


0.1228 


55 


(1,0.5,0.4) 


0.1368 


0.1982 


56 


(0.1,0.75,0.3) 


0.1354 


0.1024 


57 


(3,0.5,0.3) 


0.1342 


0.204 


continued on next column 



continued from previous column 


Config. 
index 


(e/P, T err /P, j//P)-triplet 


Distance 


Bandwidth 
savings 


58 


(2,0.5,0.3) 


0.1331 


0.2021 


59 


(0.1,0.25,0.5) 


0.1324 


0.1225 


60 


(0.5,0.25,0.3) 


0.1318 


0.1187 


61 


(1,0.5,0.3) 


0.1236 


0.1832 


62 


(0.1,0.25,0.3) 


0.1206 


0.1021 


63 


(0.1,0.75,0.1) 


0.115 





64 


(3,0.75,0.1) 


0.115 





65 


(1,0.75,0.1) 


0.1145 





66 


(2,0.75,0.1) 


0.1132 





67 


(0.5,0.75,0.1) 


0.1127 





68 


(0.5,0.5,0.8) 


0.1093 


0.1677 


69 


(0.5,0.5,0.4) 


0.1066 


0.1606 


70 


(0.5,0.5,0.5) 


0.1053 


0.1569 


71 


(1,0.5,0.2) 


0.0987 


0.1531 


72 


(2,0.5,0.2) 


0.0987 


0.1531 


73 


(0.5,0.5,0.2) 


0.0986 


0.1531 


74 


(0.1,0.5,0.2) 


0.0976 


0.1531 


75 


(3,0.5,0.2) 


0.0976 


0.1531 


76 


(0.1,0.5,0.8) 


0.0957 


0.1531 


77 


(0.1,0.5,0.4) 


0.0952 


0.1531 


78 


(0.5,0.5,0.3) 


0.0942 


0.1341 


79 


(0.1,0.5,0.5) 


0.0849 


0.1232 


80 


(0.1,0.5,0.3) 


0.0782 


0.1027 


81 


(0.1,0.25,0.1) 


0.056 





82 


(3,0.25,0.1) 


0.056 





83 


(2,0.25,0.1) 


0.0548 





84 


(0.5,0.25,0.1) 


0.0543 





85 


(1,0.25,0.1) 


0.0532 





86 


(2,0.5,0.1) 


0.0422 





87 


(0.1,0.5,0.1) 


0.0408 





88 


(3,0.5,0.1) 


0.0408 





89 


(0.5,0.5,0.1) 


0.0403 





90 


(1,0.5,0.1) 


0.0401 






TABLE I. Mapping from configuration index to parameter triplet 



